Low-loss negative refraction by laser induced magneto-electric cross-coupling 
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We discuss the feasibility of negative refraction with reduced absorption in coherently driven 
atomic media. Coherent coupling of an electric and a magnetic dipole transition by laser fields 
induces magneto-electric cross-coupling and negative refraction at dipole densities which are consid- 
erably smaller than necessary to achieve a negative permeability. At the same time the absorption 
gets minimized due to destructive quantum interference and the ratio of negative refraction index 
to absorption becomes orders of magnitude larger than in systems without coherent cross-coupling. 
The proposed scheme allows for a fine-tuning of the refractive index. We derive a generalized 
expression for the impedance of a medium with magneto-electric cross coupling and show that 
impedance matching to vacuum can easily be achieved. Finally we discuss the tensorial properties 
of the medium response and derive expressions for the dependence of the refractive index on the 
propagation direction. 

PACS numbers: 



I. INTRODUCTION 

Negative refraction of light, first predicted to occur in 
materials with simultaneous negative permittivity and 
permeability in the late 60's has become one of the 
most active fields of research in photonics in the last 
decade. Since the theoretical proposal for its realiza- 
tion in mcta-materials 0, Q and its first experimental 
demonstration Q for RF radiation, substantial techno- 
logical progress has been made towards negative refrac- 
tion for shorter and shorter wavelengths [a 111- This in- 
cludes various approaches based on split-rin g re sonator 
meta- materials 

[IHSU, photonic crystals [jl, 03, 03 
as well as more unconventional designs like double rod 
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[13, IU, Hal or fishnet structures [16 

Despite the wide variety of implementations a ma- 
jor challenge is the large loss rate of these materials 
[IU . Especially for potential applications such as sub- 
diffraction limit imaging (l9| or electromagnetic cloaking 
[20l . l2ll . I22I the suppression of absorption proofs to be 



crucial [23l. |24| . The usually adopted figure of merit 



FoM 



Re[n] 
Imfnl 



(1) 



reaches only values on the order of unity in all current 
meta-material implementations [f| with a record value of 
FcM = 3 [25[. This means that the absorption length of 
these materials is only on the order of the wavelength. 

Recently we have proposed a scheme in which coherent 
cross-coupling of an electric and a magnetic dipole reso- 
nance with the same transition frequencies in an atomic 



system [26[ leads to negative refraction with strongly sup- 
pressed absorption [27] due to quantum interference ef- 
fects similar to electromagnetically induced transparency 
[28l. l29l|. Furthermore, the value of the refractive index 
can be fine-tuned by the strength of the coherent cou- 
pling. In the present paper we provide a more detailed 
analysis of this scheme. In particular we will discuss un- 
der what conditions a magneto-electric cross coupling can 
induce negative refraction in atomic media without re- 
quiring Re[/x] < 0. The model level scheme introduced 
in [27] ] will be analyzed in detail. Explicit expressions for 
the susceptibilities and cross-coupling coefficients will be 
derived and the limits of linear response theory explored. 
The important issue of non-radiative broadenings and 
local field corrections to the response will be discussed. 
Furthermore we will derive an explicit expression for the 
impedance of a medium with magneto-electric cross cou- 
pling and show that it can be matched to vacuum via 
external laser fields such that reflection losses at inter- 
faces can be avoided. Finally we will give full account of 
the tensorial properties of the induced magneto-electric 
cross coupling and the resulting refractive index in the 
model system. We will show that in the model system 
discussed it is in general only possible to obtain isotropic 
negative refraction in 2D. 



II. FUNDAMENTAL CONCEPTS 

In the following we discuss the prospects of nega- 
tive refraction in media with magneto-electric cross- 
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coupling. The electromagnetic constitutive relations be- 
tween medium polarization P or magnetization M and 
the electromagnetic fields E and H are usually expressed 
in terms on permittivity and permeability only. Impor- 
tant aspects of linear optical systems such as optical ac- 
tivity, which describes the rotation of linear polarization 
in chiral media cannot be described in this way however. 
The most general linear relations that also include these 
effects read 



M =^E 

47T 



4tt 

" XmH. 



(2) 



Eqs.((21) describe media with magneto- electric cross cou- 
pling which are also known as bianisotropic media [3tj| . 
Here the polarization P gets an additional contribution 
induced by the magnetic field strength H and likewise the 
magnetization is coupled to the electric field component 
E. £eh and £he denote tensorial coupling coefficients 
between the electric and magnetic degrees of freedom, 
while e = 1 + 47TXe arL d p, = 1 + 47rx m arc the complex- 
valued permittivity and permeability tensors. As we use 
Gaussian units the coefficients Xe, Xm, £,eh, and £,he are 
unitlcss. 

The propagation properties of electromagnetic waves 
in such media are governed by the Helmholtz equation 



(A 



EH 



-kx)/i 1 (-k x -£he] 

UJ UJ 



E = 0. 



(3) 



A general solution of this equation for the wave vec- 
tor k is very tedious [3lj and a comprehensive discus- 
sion of the most general case almost impossible. For 
the sake of simplicity we therefore assume the permit- 
tivity e and the permeability fi to be isotropic e = el, 
fx = fxl. We furthermore restrict ourselves to a one- 
dimensional theory by choosing the wave to propagate 
in the z-direction which leaves only the upper left 2x2- 
submatrices of the tensors ^eh and £,he relevant. At 
this point we restrict the discussion to media which al- 
low for conservation of the photonic angular momen- 
tum at their interfaces. In particular, we assume the 
response matrices £eh and £,he to be diagonal in the 
basis {e + ,e_,e z } . Here e± denote circular polarization 
basis vectors e± = (e^ ± ie y )/\/2. This leads, e.g., for 
€eh, to 

£EH = £,EH e + ® e + + ^EH e ~ ® e - + €]EH e * ® e ^ 

( ^eh + Ceh)/^ -^eh-Ceh)I1 \ 
^eh-Ceh)/* {Zeh+Ceh)/* 
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(4) 



Note that the tensor @ includes also biisotropic media 
first discussed by Pendry [32| for negative refraction as 
a special case. Such biisotropic media display a chiral 



behavior, i.e., isotropic refractive indices which are dif- 
ferent for the two circular eigen-polarizations. In section 
I VII we will give an example which implements a polariza- 
tion independent but anisotropic index of refraction. 

Using ([4]) the propagation equation for a left circular 
polarized wave traveling in the z-direction can be ex- 
pressed as 



EH 



i-k: 

UJ z 



HE 



i-kl 

UJ 







(5) 



which can be solved for As fc+ is related to the 

corresponding refractive index via n + = k^c/uj we find 



n + = y e» - ttk±thL + » ^ _ £ E ) . (6) 

Note that for a vanishing cross-coupling (£eh = (,he = 
0) this simplifies to the well known expression n = y/ejl. 

Equation (O has more degrees of freedom than the 
non-chiral version and allows for negative refraction with- 
out requiring a negative permeability For example if we 
set £eh = —£,he = ££) with £ > the refractive index 
© reads 



n = yjep - i 



(7) 



Here and in the following we drop the superscript + for 
notational simplicity. 

Magnetic dipole transitions in atomic systems are a 
relativistic effect and thus magnetic dipoles are typically 
smaller than electric ones by factor given by the fine 
structure constant a w 1/137. Since furthermore mag- 
netic resonances are typically not radiatively broadened 
the magnetic susceptibility Xm per dipole is typically less 
than the electric susceptibility Xe by a factor given by the 
fine structure constant squared a 2 . On the other hand, as 
we will show later on, the cross-coupling coefficients £eh 
and £he scale only with one factor of a. Thus Eq. (|7|) 
represents an improvement compared to non-chiral ap- 
proaches, since a negative index can be achieved at den- 
sities where [i is still positive but £ > ^/ejl. Negative 
refraction with Re[/i] > comes with the requirements 
to find large enough chiral coupling coefficients ^eh and 
£,he and additionally to control their phases in order to 
get close enough to £eh = —£,he = i£ necessary for neg- 
ative refraction. 

As in [23] we will analyze these fundamental concepts 
in more detail and consider a modified V-type three-level 
system (Fig. [I]). It consists of an electric dipole transi- 
tion which couples the ground state |1) and the excited 
state 1 3) as well as a magnetic dipole transition between 
|1) and the upper state |2). We assume states |2) and |3) 
to be energetically degenerate such that the electric (E) 
and magnetic (B) field components of the probe field can 
couple efficiently to the transitions |1) — |3) and |1) — |2), 
respectively. In order to couple the electric and mag- 
netic dipoles, i.e., to induce a cross-coupling in the sense 
of equation we add a strong resonant coherent cou- 
pling field between the two upper states |2) and |3) with a 
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FIG. 1: A simplistic three-level system which employs electro- 
magnetically-induced cross-coupling. E, B are the electric 
and magnetic components of the probe field, 72 and 73 the 
decay rates out of levels |2) and |3). Q c is an applied field 
which couples levels |2) and |3). 




FIG. 2: (color online) Modification of the level scheme of 
Fig. [T] The ground state |1) is substituted by the dark state 
\D) = (fi 2 |l) - nil^/^/fif+nf of the 3-level A-type sub- 
system formed by levels {|1}, |4), |5)}. 



Rabi-frequency fi c . Note that the configuration of Fig. Q] 
complies with the requirements of parity rules. The mag- 
netization of the system at the probe field frequency is 
given by the coherence of the |1) — |2) transition which 
also gets a contribution induced by the electric field E 
and likewise the polarization, connected to |1) — |3), is 
not only induced by E but also by the magnetic field B. 
Following the discussion above the level scheme of Fig. [1] 
should therefore show a negative refractive index without 
requiring Re[ju] < 0. 

Concerning the absorption we note that the main con- 
tribution to the imaginary part of n stems from the per- 
mittivity e. The radiative decay rate 72 of the magnetic 
dipolc transition |2) — |1) is typically a factor a 2 smaller 
than 73, the decay rate of the electric dipole transition. 
As a consequence the strong field fi c couples state |3) 
strongly to the metastable state |2) which, on two-photon 
resonance, is the condition for destructive quantum inter- 
ference for the imaginary part of the permittivity, known 
as electromagnetically induced transparency (EIT) [28|. 
Additionally, for closed loop schemes (Fig.QJ it is known 
from resonant nonlinear optics based on EIT [H, [34|, 
that the dispersive cross-coupling, which in our case is the 
magneto-electric cross coupling, experiences constructive 
interference. 

In essence the coupling of an electric to a magnetic 
dipole transition should lead to negative refraction for 
significantly smaller densities of scatterers compared to 
non-chiral proposals [H, [H, ■ We additionally expect 
the imaginary part of e, which represents the major con- 
tribution to absorption, to be strongly suppressed due 
to quantum interference effects while simultaneously the 
cross-coupling terms should be further enhanced. 

Though conceptually easy the scheme of Fig.Q]has sev- 
eral drawbacks which demand a modification of the level 
structure, (i) As stated above the phase of the chiral- 
ity coefficients £eh, £,he must be adjustable to control 
the sign of the refractive index and induce Re[n] < 0. 
As fl c is a dc-field in the scheme of Fig. [1] the phases of 



£eh, £,he are solely given by the intrinsic phase of the 
transition moments and therefore can not be controlled, 
(ii) To suppress the absorption efficiently there must be 
high-contrast EIT for the probe field. The critical pa- 
rameter for this effect is the dephasing rate 721 of the 
coherence p 2 i between the two EIT "ground" -states |2) 
and |1). Since state |2) has an energy difference to state 
|1) on the order of the probe field frequency, the coher- 
ence P21 is highly susceptible to additional homogeneous 
or inhomogeneous broadenings which ultimately can de- 
stroy EIT. (iii) The level structure must be appropriate 
for media of interest (atoms, molecules, excitons, etc.). 
Although the scheme of Fig. Q] is not forbidden on fun- 
damental grounds it is very restricting to require that 
electric and magnetic transitions be energetically degen- 
erate while having a common ground state. 

One possible alternative level structure which solves 
these problems is shown in Fig. [5] The former ground 
state |1) is now substituted by the dark state \D) of 
the 3-level A-type subsystem formed by the new lev- 
els {|1), |4), |5)} as shown in Fig. \D) is determined 
by the two coupling fiel d Rabi fr equencies and f^: 
\D) = (fi 2 |l) - Oi|4))/i/Of+T2f. This simple manipu- 
lation indeed solves the above mentioned problems: The 
upper states |2) and |3) are no longer degenerate, i.e. 
the coupling Rabi frequency £l c is now given by an ac- 
field. By adjustment of f2 c the phase of ^eh and £he 
can be controlled. The critical parameter of EIT for this 
scheme is the dephasing rate 724 of the coherence p2i 
between states |2) and |4). By assuming levels |2) and 
|4) of Fig. [2] to be close to degenerate p 2 4 can be taken 
to be insensitive to additional broadenings. Compared 
to the 3-level system of Fig. [T] the electric and magnetic 
transitions here do not share a common state while the 
transition frequencies are still degenerate. This leaves 
much more freedom regarding a realization in real sys- 
tems. 

In the next section we will analyze the 5-level-system 
in detail. 
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FIG. 3: 5-level scheme for the implementation of negative re- 
fraction via electromagnetically induced cross-coupling. The 
magnetic dipole transition |2) — |1) and the electric dipole 
transition |3) — |4) are coupled by O c to induce chirality. 
The "ground" -state of the syst em is fo rmed by the dark 
state \D) = (fi 2 |l) - Sli\4))/y/Q,l + of the subsystem 
{|1),|4),|5)}. 



III. 5-LEVEL SCHEME 
A. Linear response 

The level scheme in question is shown in greater detail 
in Fig. [3] The transition |2) — |1) is magnetic dipole 
by nature, all other transitions are electric dipole ones. 
Note that this complies with the demands of parity. It 
is assumed that for reasons of selection rules or absence 
of resonance no other transitions than the ones sketched 
in Fig. [3] are allowed. The Hamiltonian H = H at om — d • 
E(t) — fi ■ B(t) of the system can be written explicitly as 



- X -d^Ee-^ |3) (4| - i^Be-*"** |2) (1| 
|0 ie -^*|5) (l|-^n 2e -^*|5) (4| 



--n c e-^*|3) (2| + h.c. 



(8) 



Here CZ34 = (3| e r ■ |4) and fi 2 i = (2| /x • ee |1) are the 
electric and magnetic dipole moments, E and B the elec- 
tric and magnetic components of the weak probe field 
which oscillates at a frequency u> p . The Rabi frequencies 
fii, and f2 c belong to strong coupling lasers which 
oscillate at frequencies u)\, L02 and u> c , respectively. We 
choose c?34 and /x 2 i as well as Oi and Q2 to be real whereas 
the strong coupling Rabi frequency Sl c has to stay com- 
plex for the closed loop scheme. 



To include losses in our description we solve for the 
steady state solutions of the Liouville equation of the 
density matrix. In doing so we introduce population 
decay rates 74, i e {1,2,3,4,5}. As we focus on the 
linear response we treat the probe field amplitudes E 
and B as weak fields which allows to neglect the upper 
state populations P33 and p22- In contrast the subsystem 
{|1), |4), |5)} contains strong fields with Rabi frequencies 
f2i and r?2 and should be treated non-perturbatively. 

We first solve the 3-level subsystem {|1), |4), |5)} undis- 
turbed by the probe field. Under the assumption that |4) 
is meta-stable and therefore 74 -C 75 holds, the exact so- 
lution is: 



P11 



~(0) 

Pal 



(0) 

P44 



Nil 



Nil 



(0) 

P55 



-(0) 

P&i 



-(0) 
= P 5 4 



(9) 



0. 



This solution for the A-type subsystem indee d corrc- 
sponds to the pure state \D) = (f2 2 |l)-Oi|4))/ v /0 1 r + fi| 
via p suhs y s - = \D) (D\. Note that density matrix compo- 
nents with tildas denote slowly varying quantities. 

We proceed by solving for the polarizabilities of the 
complete 5-level-schemc up to first order in the probe 
field amplitudes E and B. Since the induced Polarization 
P is proportional to the coherence of the electric dipole 
transition ^34 whereas the induced Magnetization M is 
proportional to the density matrix element p~2i we arrive 
at 



P =gd 34 p 34 := ga^^E + ga^^B 
M =QH2ip2i ~ ga BE E + ga BB B. 



(10) 



Here g is the number density of scatterers, a EE , a EB , 
a BE and a BB are the direct and the cross-coupling po- 
larizabilities. They are given by: 



Q 



EE 



^34/°44 (742 + i{A E - S c )) 



, BB 



as well as 



2h (742 + i{A E - <5 C ))( 7 34 + iA E ) + |ft c | 2 /4' 

(11) 

J mIiPi°i ) (73i + i(A B + 5g)) 

2h (731 + i(A B + <5 C ))( 721 + »A B ) + |« c | 2 /4 

(12) 



. EB 



.BE 



1 d 3 4p,2ipfi^c 

Ah (742 + i{A E - 5 C ))(734 + iA E ) + |ft c | 2 /4' 

(13) 



1 



<i34A*2lP41 ) ^ i * 



4S (731 + i(A B + 5 c ))(72i + iA B ) + N c | 2 /4' 

(14) 

Here the definitions 7^ = (73 + 7j)/2, A E = uj 34 — lo p , 
A B = W21 — and S c = uj 3 2 — w c apply with oj^v = 
lu^ — lu^ being the transition frequencies between levels 
\fx) and \v). Note that these solutions are only valid 
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FIG. 4: (color online) Real (solid lines) and imaginary (dashed 
lines) parts of the electric (a EE ) and magnetic (a BB ) polariz- 
abilities as well as the cross-coupling parameters (a EB , a BE ) 
for arbitrary but the same units for Q. c = 10 4 72e I7r ' /2 . 



for A E = A B 
condition 



5 C which corresponds to the resonance 



UJ C = LUl — U>2 



(15) 



which ensures the total frequency in the closed loop 
scheme to sum up to zero. 

In order to visualize the polarizabilitics we set the 
magnetic dipole decay rate to a typical radiative value 
for optical frequencies [38j of 72 = 1kHz, the electric 
dipole decay rates correspondingly to 73 = 75 — 137 2 72 
and 71 = 74 = for the (meta-)stable states |1) and 
1 4). The electric and magnetic dipole matrix elements 
(I34 and /X21 are determined from the radiative decay 
rates via the Wigncr-Wcisskopf result [3!| ^34(^21) = 
\Z373(72)fi.c 3 /(4kj 3 ) for an optical frequency correspond- 
ing to A = 600 nm. The Rabi frequencies of the A- 
type subsystem attain real values f2i = Q2 = 10 2 72 
while the coupling Rabi frequency Q c is chosen com- 
plex, f] c = |n c |e z ^. We furthermore specialize to S c = 
which implies = A^. In order to have increas- 
ing photon energy from left to right in figures presented 
here all following spectra are plotted as a function of 
A = -A E = -A B . 

For f2 c = the cross-coupling coefficients a EB and 
a BE vanish exactly whereas the electric as well as the 
magnetic polarizability both show a simple Lorentzian 
resonance. Introducing a non-zero coupling strength fi c 
changes the response dramatically. As shown in Fig. 0] 
for fl c = 10 4 72e l7r / 2 the electric polarizability a EE in- 
deed displays clectromagnctically induced transparency 
(EIT). The value of \fl c \ is optimized as will be dis- 
cussed in section ITlI CI As long as the coupling field fi c is 
present Im^^^] on resonance is proportional to the de- 
coherence rate 742 of the two EIT ground states |4) and 
1 2) which can become very small. Thus the prominent 
feature of EIT emerges: suppression of absorption. In 
contrast the magnetic polarizability a BB still shows an 
ordinary Lorentzian resonance since 731 is always large. 
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FIG. 5: (color online) Deviation (|16[1 of real and imaginary 
parts of the exact polarizabilities a IJ (E, B) compared to the 
linear response results for 2 different probe field Rabi frequen- 
cies: Qe = 137 Sis = 72 (real: solid lines, imaginary: dashed 
lines) and Qe = 137 Sis = 10-72 (real: dash-dotted lines, 
imaginary: dotted lines). 



Due to the coupling to the strong electric dipole transi- 
tion the magnetic resonance is broadened as compared 
to its radiative lincwidth which is accompanied by a sig- 
nificant decrease of the magnetic response susceptibility 
on resonance. 

For a non-vanishing $l c the two cross-couplings a EB 
and a BE show strongly peaked spectra. Note that the 
phase 4> = 7r/2 has been chosen such that on resonance 
a EB = —a BE oc i holds, as demanded in section ITT1 

Note furthermore that we verified numerically that all 
polarizabilities and cross-coupling terms are causal and 
thus fulfill Kramers-Kronig relations. 



B. Limits of linear response theory 

For radiatively broadened systems 72^ ~ « 2 72i^ 
holds. Thus magnetic transitions saturate at much lower 
field amplitudes than corresponding electric dipole tran- 
sitions. For this reason we have to analyze the saturation 
behavior of the system. 

To rule out saturation effects and validate the use of 
linear response theory we solve the Liouvillian equation 
for the 5-level system of Fig. [3] to all orders in the electric 
and magnetic field amplitudes E and B which can only 
be done numerically. We determine the polarizabilities 
a EE (E,B), a BE (E,B), a BB {E,B) and a EB {E,B) from 
the numerically accessible density matrix elements P34 
and P21 (see |Appendix A[ ) for the worst case scenario of 
purely radiatively broadened transitions and compare to 
the results of the linear response theory. Fig. [5] shows 
the deviation of the exact solution to the linear response 
results 
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(16) 
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for all polarizabilitics. The solid lines correspond to 
£ = 137 H B = 72, the dashed lines to fl E = 137 = 
10 -72- Here He = d^E/h and Hb = ^lB/h denote the 
electric and magnetic probe field Rabi frequencies for the 
5-level system. 

For a radiatively broadened electric dipole 2-level atom 
we can estimate the probe-field Rabi frequency required 
to lead to a 1% upper state population to be He — 
VO.OI73 f» 13.772. ^From Fig. [5] we note that even for 
He = 137 0b = IOX72, i.e., the same order of magnitude 
as He, the deviation of the exact result from the spec- 
trum obtained in linear response approximation never 
exceeds 10 _2 . As a result we conclude that the 5-level 
scheme is not significantly more sensitive to saturation 
effects than any ordinary electric dipole transition. 

This behavior is a result of the coupling Rabi frequency 
H c due to which the \2) — |1) transition experiences an 
additional broadening which makes it less susceptible to 
saturation. 
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FIG. 6: (color online) Real (solid lines) and imaginary (dashed 
lines) parts of the electric (a EE ) and magnetic (a BB ) polariz- 
abilities as well as the cross-coupling parameters (a EB , a BE ) 
for arbitrary but the same units for Q c = 10 4 7 2 e i7r/2 . In con- 
trast to Fig. [4] additional homogeneous broadenings according 
to Eq. IQU with 7 P = 10 3 7 2 apply. 



C. Non-radiative broadenings 

As noted in section [J additional broadenings are essen- 
tial for the description of spectral properties as soon as 
magnetic transitions are involved. 

To add an additional homogeneous dephasing rate we 
formally add to every static detuning A^, As, S c an 
extra random term xe, xb, x c (e.g. A^ — ► A^ + xe) 
which has to be convoluted with a Lorcntzian distribution 



L{x) = - 



1p 



1 



7T J 2 + X 2 ' 



(17) 



For, e.g., a EE this amounts to 



a EE 



dxE<ixBdx c a EE (xe, xb, x c )L(xe)L(xb)L(x c 



(18) 

Levels |2) and |4) are approximately degenerate and thus 
assumed to experience correlated phase fluctuations. As 
a consequence 742, which is relevant for EIT, remains 
unchanged while the same width 7 P applies for both A^ 
and S c . For simplicity we choose the same width for A^. 
The convolution integral can be solved analytically which 
results in the substitution of the off-diagonal decay rates 
Hj, i ^ j (fTT|) - (fT4l) according to 



|fi c | = 10 4 72- Since 742 remains unbroadened the elec- 
tric polarizability a EE still shows EIT while the spec- 
trum of a BB shows a simple but broadened resonance 
with a EE w 137 2 a BB . Similarly a EE w 137a BB and 
a EE w l37a BE hold approximately 

To incorporate the effect of an inhomogeneous Dopplcr 
broadening mechanism on the spectrum the same formal- 
ism as for the homogeneous case, but with a Gaussian 
instead of a Lorcntzian distribution, can be used. 



D. Local field correction 

So far we have dealt with a single individual radia- 
tor (e.g., atom) responding to the locally acting fields. 
For large responses these local fields are known to differ 
from the averaged Maxwell field. We therefore have to 
correct the results of section UlI Al bv the use of Clausius- 
Mossotti type local field corrections. Note that due to 
the cross-coupling the influence of the magnetic proper- 
ties is enhanced by a factor of approximately 137. We 
therefore also include magnetic local field corrections in 
the treatment. 

To add local field corrections to the response the fields 
E and B Eq. (fT0|) are interpreted as local or microscopic 



742 
734 



742 
734 



721 
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■27 P - 



(19) 



In contrast 721 and 734 encounter a broadening j p , like- 
wise 731 experiences a broadening 27 p . 

We choose the value 7 P = 10 3 72 which is typical for 
rare-earth doped crystals at cryogenic temperatures. For 
a given broadening "f p the cross-coupling terms reach a 
maximum for the coupling Rabi frequency attaining the 
optimal values 13772 and 1370v / 2072> respectively. Fig. [(J 
shows the polarizabilities for an intermediate value of 
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P =g a BB £ mlcro + ga^B" 
M =ga BE E micm + ga BB B micl ° 



(20) 



(g: number density of scatterers). The relations between 
the local and the corresponding macroscopic field ampli- 
tudes can be obtained from phcnomenological considera- 
tions WO, WlJ which read 
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FIG. 7: (color online) Real (solid lines) and imaginary (dashed 
lines) parts of the refractive index, including local field effects 
for two different densities. 



for the electric and the magnetic field, respectively. Note 
that we need to replace B by H to find the permittiv- 
ity e, permeability fi, and coefficients £eh and £he of 
equation ([2|) in terms of the polarizabilities of equations 
(fTTj) - (fT4|). This can be done most easily for the local 
microscopic fields for which B micl '° = H micro holds [Hj]. 
Solving Eq. ([2(7)) together with (|2"Tj) for P and M in terms 
of the macroscopic field amplitudes E and H yields 



e =1 + 47T— t— x 

£loc 



_ EB „BE EE „BB 

g\ a a —a a 



u=l + An— ^— x 



(22) 



ten =^^ EB , 
^E=^ c a BE , 
with the denominator 

C^=l-^ga EE -^ga BB 
'4tt n 



(23) 



„EB„BE EE „BB 

a a — a a 



)■ 



Note that for media without a magneto-electric cross- 
coupling a rigorous microscopic derivation [42j validates 
the phenomenological procedure adopted here. 



E. Negative refraction with low absorption 

With the permittivity e and the permeability /i given 
by Eq. ([22]) and the parameters £eh and £he [Eq. ((23)) ] 
we determine the index of refraction from Eq. (|6|) . As an 
example, Fig. [TJa) shows the calculated real and imag- 
inary parts of the refractive index as a function of the 
probe field detuning A for a density of g = 5 • 10 14 cm -3 
and otherwise using the parameter values defined in sec- 
tions IIII Al and IIII CI The shape of the spectrum is gov- 
erned by the permittivity e with the prominent features 



FIG. 8: (color online) (a) Real (solid line) and imaginary 
(dashed line) parts of the refractive index for g = 5T0 16 cm~' ! . 
Compared to Fig. [7jb) here = applies, (b) Real (solid 
line) and imaginary (dashed line) parts of the refractive index 
as a function of the phase <f> of the coupling Rabi frequency 



of EIT: suppression of absorption and steep slope of the 
dispersion on resonance. Clearly for this density there is 
no negative refraction yet. 

In Fig. [7Jb) the spectrum of n for an increased den- 
sity of g = 5 • 10 16 cm~ 3 is shown. Note that in contrast 
to Fig. [TJa) the frequency axis is scaled in units of the 
broadening -f p as local field effects start to influence the 
shape of the spectral line at this density. We find sub- 
stantial negative refraction and minimal absorption for 
this density. The density is about a factor 10 2 smaller 
than the density needed without taking chirality into ac- 
count f35l ]. 

In order to validate that the negative index results from 
the cross-coupling we compare the spectrum of Fig. [TJb) 
to a non-chiral version. As setting Sl c = influences the 
permittivity and permeability as well, we set = 



by hand such that the cross-couplings vanish identically. 
The resulting index of refraction is shown in Fig. [SJa) 
for a density g = 5 ■ 10 16 cm~ 3 . We find that without 
cross-coupling no negative refraction occurs. Thus the 
negative index at this density is clearly a consequence of 
the cross-coupling. 

Following the qualitative discussion in section |TT] we 
have set the phase <j> of the coupling Rabi frequency tt c to 
<f) = tt/2. Fig.[8jb) shows the dependence of the refractive 
index on <fi taken at the spectral position A = — 0, 0457 P 
at which Re[n] reaches its minimum for g = 5 ■ 10 16 cm~ 3 . 
As expected the refractive index is strongly phase de- 
pendent. A change of the phase by 5(f) = it for example 
reverses the influence of the cross-coupling and gives a 
positive index of refraction Re[n] > 0. Note that the sym- 
metry Re[n((/>)] = — Re[n(27r — <fi)] is coincidental since for 
the chosen parameters e ~ 0. 

By increasing the density of scatterers g further the 
optical response of the medium increases. For a spec- 
tral position slightly below resonance (A = —0, 0457 P ) 
where negative refraction is obtained most effectively we 
show Re[e] as well as Imf^Ejj] and Im[£jj£:] as functions 
of the density g [Fig. [5{a)]. Due to local field corrections 
the permittivity is of the same order of magnitude as the 
cross-coupling terms. The imaginary parts of the param- 
eters £,eh and ^he increase strongly with opposite signs 
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FIG. 9: (color online) Real part of (a) the permittivity e 
(solid line) as well as the imaginary parts of ^eh (dashed line) 
and t^HE (dotted line) and (b) real (solid line) and imaginary 
(dashed line, x 50) parts of the refractive index, as well as the 
real part of the permeability (dotted line) as a function of the 
logarithm of the density log g. 
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FIG. 11: (color online) The Figure of Merit FoM for densities 
£i = 5- 10 17 cm~ 3 (solid line) and g = 5 ■ 10 16 cm~ 3 (dashed 
line) as a function of the detuning A. 
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FIG. 10: (color online) Real (solid line) and imaginary 
(dashed line) parts of the refractive index, including local field 
effects for a density of g = 5 • 10 17 cm -3 . 



causing the refractive index to become negative. The 
corresponding density dependence of the refractive index 
is shown in Fig. H^b). We find that Re[n] becomes nega- 
tive while the absorption Im[n] stays small (note that in 
Fig. E^b) lm[n] is amplified by a factor of 50). Addition- 
ally Rc[/x] is positive and becomes larger for increasing 
density as a consequence of operating on the red detuned 
side of the resonance (A < 0). 

As an example for higher densities the spectrum of n is 
shown for g = 5 ■ 10 17 cm~ 3 in Fig. [TO] Compared to the 
case of q = 5 • 10 16 cm~ 3 [Fig. HJb)] Im[n] did not change 
much qualitatively while Re[rt] reaches larger negative 
values. 

Remarkably, in Fig. [9jb) we find that the absorption 
Im[n] reaches a maximum and then decreases with in- 
creasing density of scatterers. This peculiar behavior is 
due to local field effects which invariably get important 
at such high values of the response. As a consequence the 
spectral band with minimal absorption broadens with in- 
creasing density due to local field effects. Hence the cho- 
sen spectral position moves from the tail of the band edge 
to the middle of the minimal absorption band. 

As a consequence of the low absorption and corre- 
sponding increasing values of Re[n] the FoM continues 
to increase with density and reaches rather large val- 



ues. In Fig. [TT] we show the FoM as a function of A 
for g = 5 ■ 10 16 cuT 3 and g = 5 ■ 10 17 cm" 3 . While 
the FoM reaches for g = 5 • 10 16 cm~ 3 values of as 35 
it climbs for g = 5 • 10 17 cm~ 3 up to FoM ps 350. These 
results should be contrasted to previous theoretical pro- 
posals and experiments on negative refraction in the op- 
tical regime, for which the refraction/absorption ratio is 
typically on the order of unity. 

IV. TUN ABILITY 

As first noted by Smith et al. [H| and Merlin pi ] 
sub-wavelength imaging using a flat lens of thickness d 
requires not only some negative refractive index but an 
extreme control of the absolute value of n. For an in- 
tended resolution Aa; the accuracy with which the value 
n = — 1 (assuming a vacuum environment) has to be met 
is given by 

|An|=exp|-^}. (24) 

For a metamaterial with Re[n] < Eq. (|24|) presents 
a considerable obstacle for the operation of a superlens 
approaching far field distances (d 3> A) as it demands 
an extreme fine-tuning of the refractive index in order to 
achieve a resolution beyond the diffraction limit. 

Our scheme allows to achieve such a fine tuning. In 
Fig. [T^] we show the real and imaginary parts of n as 
a function of log[|fi c |/7 3 ] for a density of g = 1,56 • 
10 17 cm -3 . As the coupling Rabi frequency approaches 
73 we find small Im[n] while Re[n] attains negative val- 
ues. The dispersion then changes only slightly with a 
small slope and values around n = — 1. Therefore the 
refractive index can be fine tuned by relatively coarse ad- 
justments of the strength of the coupling field f2 c . Note 
that the value and slope of n for |fi c | f=a 73 can be cho- 
sen by adjusting the density g of the medium and the 
spectral position of the probe light frequency. 

Apart from potential imaging applications the 5-level 
quantum interference scheme allows for devices operating 
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material equations ([2]). For medium 1 we get 

k r x e,E r e~ iklZ = 



l^xe+E^ 



iki z 



-Mi (-Hie** 1 
c 



e+. 



(28) 



Noting that e z x e± = ±ie± holds, this simplifies for 
z = to the scalar equation 



ifci (i?; — E r 



-m(Hi + H r ) 



(29) 



FIG. 12: (color online) Real (solid line) and imaginary 
(dashed line) parts of the refractive index as a function of 
the coupling field Rabi-frequency |fl c | relative to the radia- 
tive decay rate 73, for g = 1, 56 • 10 17 cm -3 . 



where k, = — k r = k\e z has been applied. Similarly we 
obtain for medium 2 



ik 2 E t = -^ HE E t + LL 2 H t ). (30) 

c 



in a wide range of positive and negative refractive indices 
with simultaneously small losses. 



V. IMPEDANCE MATCHING 

When considering the applicability of optical devices 
reflection at boundaries between different media play an 
important role. Impedance matching at these bound- 
aries is often essential for the performance of the device. 
For sub- wavelength resolution imaging the elimination of 
reflection losses is particularly important, see Eq. (|2~i)) . 
In this section we thus derive conditions under which the 
boundary between non-chiral and chiral, negative refract- 
ing media is non or little reflecting. 

We consider a boundary between a non-chiral medium 
1 (z < 0) with £1, fi\ and medium 2 (z > 0) which 
employs a cross-coupling (£2, M2, £,eh, Che)- We assume 
again a wave propagating in z-direction in a medium cor- 
responding to the tensor structure (0| such that we can 
restrict to an effectively scalar theory for, e.g., left circu- 
lar polarization. 

We decompose the e + field component in region 1 into 
an incoming Ei and a reflected part E r 



Ei(r) = (Eid 



ik±z 



E r e 



—ik\z 



(25) 



(ki = |kj| = |k r |). In medium 2 (z > 0) only a trans- 
mitted wave E t exists due to the boundary condition at 
infinity 



E 2 (r) 



E t e ik2Z i 



(26) 



The connection of these modes at the interfaces and 
similar ones for the magnetic field H(r) is established 
by the boundary conditions n x (E2 — Ei) = and 
n x (H 2 - Hi) =0. At z = we find 



Ei + E r 



E t 



Hi + H r = H t 



(27) 



Moreover an independent set of equations is found from 
Maxwell's equations in Fourier space together with the 



We eliminate the magnetic field amplitudes from (|27|) . 
(|29|) . and ([30]) and solve for the ratio of reflected and 
incoming electric field amplitudes which yields 



Er 

Ei 



Mi n 2 + i£,HE 

M2 



Mi n 2 + i£he 

M2 




(31) 



Here the wave numbers k\ and k 2 have been replaced by 
k\ = niLu/c = ^jE\[i\ (jj/c and k 2 = n 2 u>/c. Equation ((3"Tj) 
is a generalization of the well-known Fresnel formulas for 
normal incidence to a cross-coupled medium. Impedance 
matching is defined as the vanishing of the reflected wave 
E r = 0, i.e., a complete transfer of the incoming field into 
medium 2: 



fjiin 2 + i£ H E 

M2 



(32) 



Using the explicit form of n 2 for the particular polariza- 
tion mode © we find the more convenient expression 



(,EH + S.HE 
2^2M2 



1 £,EH + £HE 

2 ^£2^2 



(33) 

which obviously simplifies to the well known result 
\/ei/fii = \fe 2 f[i 2 for the non-chiral case £eh = Che = 
0. The right hand side of (|33[) is the inverse impedance 



of the cross-coupled medium. For causality reasons 
the square root in Z 2 for passive media has to be taken 
such that Re[Z2] > is obtained [43|. Fig. [TBI shows the 
real and imaginary parts of Z^ 1 and of the index of re- 
fraction. For the case s\ = /xi = 1 the impedances of 
the two media are matched at the interface as soon as 
Z^ 1 = 1 + iO holds. Applying the density g = 1.56 • 10 17 
3 we find Z^ 1 



cm 



1.003 

-2 



z0.0006 at the spectral po- 
sition A = 1.17 ■ 10~ 2 -f p while the index of refraction 
attains n = —1.0003 + i0.009. The corresponding figure 
of merit is about FoM m 110. 
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FIG. 13: (color online) Real (solid line) and imaginary 
(dashed line) parts of the refractive index as well as the real 
(dash-dotted line) and imaginary (dotted line) parts of the 
impedance Z^ 1 for g = 1, 56 • 10 17 cm -3 . 
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FIG. 14: (color online) For a direction of incidence other than 
the direction of the coupling field vector E c numerous addi- 
tional angle dependent couplings occur. 



VI. BEYOND ID: ANGLE DEPENDENCE 



In section [TTI we specialized our discussion to an effec- 
tively scalar theory by restricting to a particular direc- 
tion of propagation and left circular polarization (e+). 
We now want to analyze the dependence of the refractive 
index on the propagation direction of the light, which 
requires to take into account the tensor properties of all 
linear response coefficients. To this end we consider a 
generalization of the 5-level scheme Fig. [3] that includes 
the full Zeeman sublevel structure shown in Fig. [TU The 
coupling field E c is assumed to be linear polarized in the 
z-direction. As the quantization axis we choose the prop- 
agation direction of the probe light. In this scheme the 
requirements of section [TT] are fulfilled. 

Let us first consider the case when the probe light prop- 
agates along the z axis. Due to Clebsch-Gordon rules 
this solely leads to couplings and between the 
transitions |3, +) — |2, +) and |3, — ) — |2, — ), respectively. 
As a result the two circular polarizations of a probe field 
traveling in z-dircction are eigenmodes and the scalar 
treatment from section [TTI is valid. We therefore get 



section UlI CI we find the relations 



4 ± 2 (*eh - She) ( 34 ) 



for the left [cf. Eq. ((6]))] and right circular polarizations. 
The Wigner-Eckart theorem [38[ implies that the electric 
dipole moments of the |3, +} — |4) and the |3, — ) — |4) 

Similarly the matrix ele- 



transition coincide: 



"34 



*34- 



mcnts of the magnetic dipole transitions are independent 
of the polarization state: fi^i = Mil- ^ n contrast we find 



"32 



d 32 . Thus the coupling Rabi frequencies of the 

^ 'E c |/ft have 



left and right circular branches fij^ = d 32 
a relative sign 



o++ = (35) 

From (jlip - (fl~4]) together with the results obtained in 



=e 



11' = fl 

£~HE = ~^HE 



(36) 



One recognizes that the refractive indices of e + - and e_ 
polarizations arc identical 



(37) 



Hence the index of refraction of the scheme from Fig.fTHis 
independent of the polarization state of probe light prop- 
agating in z-direction. For this reason the electromagnet- 
ically induced cross-coupling in the scheme of Fig.fMldoes 
not correspond to a chiral medium for which the circular 
components should have different refractive indices. 

Next we allow for an angle between the propagation 
direction of the probe field and the direction of the cou- 
pling field vector E c . In particular we employ a frame 
of reference whose z-direction is fixed to the k-vector: 
k ~ e z . The direction of the coupling field which is fixed 
with regard to the laboratory frame is then given by polar 
angles 6, <j>: 

(E CtX ,E c ,y,E c . z ) = |E c |(sin0cos0,sin0sin(/>,cos0). 

(38) 

As the atomic quantization axis is assumed to be given by 
the e z -axis of the k-framc, the probe field will encounter 
an unchanged atomic level structure irrespective of the 
direction of propagation. As indicated in Fig. [M] the 
|1) — 1 4} transition is assumed to be a J = 0,M = 
to J = 0, M = transition and thus the dark state is 
spherically symmetric and does not depend on the polar 
angles and (j>. 

In this framework angle-dependent propagation is tak- 
en into account by means of angle dependent coupling 
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Rabi frequencies 



n 



-o 



f2 Ci o cos 6 

71" 
7f 



= -n~, 

sin<9e~^ = 



(39) 
(40) 

(41) 



according to (J38|). Here fl cfl = (3||d||2) • |E c |/(\/6S) is 
found from the Wigner-Eckart theorem where (3||d||2) 
denotes the reduced dipole matrix element. 

The angle-dependent cross coupling tensor for the elec- 
trically induced magnetization reads 



3 IT 



FIG. 15: (color online) Real (solid line) and imaginary 
(dashed line) parts of the refractive index (|44[) as a function 

of e. 



( 



a BE = a BE 



cos 9 




sin( 



cos v sin ( 



-ic/> 



V2 

e +i<j> 



e +i<i> e" 
smf) — =r- sinf 



\ V2 V2 



(42) 



/ 



with a BE given by (fTl| . Note that the tensor is expressed 
in the {+, — , z}-basis. For example the coefficient a EB 
which describes the e + -polarized electric field induced by 
a e z -polarized magnetic field (in the k-frame) is given by 
the upper right entry. For a EB we find (|42[) as well but 
with a BE replaced by a EB from equation (fl~3|) . On the 
other hand the electric polarizability is given by 



a EE = a EE l 



sin 



a EE \n c \ 2 

D42D34 



sin 2 9 e 2 ^ 



sin 2 9 e" 2 ^ 



sin0cos0e-^ sin 20 



sin ft cos ft e 

sin2fle-^ 
8\/2 
cos 2 



V 4V2 



8^2 



(43) 



with a EE determined by Eq. (jTTJ) and D42 = (742 + 
i(As — 5c)), -D34 = (734 + «Ab). For the magnetic polar- 
izability a BB the same tensor structure applies. Again 



a has to be replaced by (fT2jl and -D42-D34 is substituted 
by D31-D21 = (731 + + 5 c ))(72i + a'A B ). 

For incidence in the z-dircction the tensors simplify 
significantly. The cross-couplings reduce to a tensor pro- 
portional to e + <£) e*^_ — e <E) e*_ which identically corre- 
sponds to (J4j) for £ z = and £ + = — In the same limit 
(9 = 0) the electric and magnetic polarizabilities become 
diagonal. In particular a EE and a EE are given by (fTTj) 
and therefore potentially display EIT while the a EE en- 
try simplifies to a simple Lorentzian resonance structure. 
In contrast the diagonal elements of a BB always display 



a Lorentzian resonance with (a 
(a BB ) coherent coupling. 



BB 

+ ' 



a__) and without 



From the angle dependent response tensors we find an 
angle dependent index of refraction. The true index of 
refraction which takes the full form of (|4"3"|) into account, 
i.e., the angle dependent correction to e and [A, gets very 
complicated. We here note that under the assumption of 
isotropic permittivity and permeability we find the fairly 
simple result 



EH 



- £,eh£,he 



+ t;(£eh - £he)cos(9). 



HE 



] COS 2 (9) 



(44) 



independent of the polarization state. We conclude that 
even the idealized case e ~ 1 ~ p, does not give an 
isotropic index of refraction. 

In Fig. [15] we show the index of refraction (|44j) clS cX 
function of the polar angle 9. We use values of the 
response functions taken at a spectral position A = 
— 0,0357 P for a density g = 5 • 10 16 cm -3 . We emphasize 
that the angle dependence in (|44j) results in an index of 
refraction which varies over a broad spectrum of positive 
and negative values for different angles. 



VII. CONCLUSION 

In conclusion we have shown that coherent magneto- 
electric cross-coupling improves the prospects to obtain 
low-loss negative refraction in several ways. The densi- 
ties needed to get Rc[n] < arc small enough to consider 
implementations in, e.g., doped crystals. The presence 
of quantum interference effects similar to electromagnet- 
ically induced transparency suppresses absorption and 
at the same time enhances the magneto-electric cross- 
coupling. As a result our scheme allows for a tunable low- 
loss negative refraction which can be impedance matched 
by means of external laser fields. 
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Appendix A. EXACT NUMERICAL SOLUTION 
OF THE LIOUVILLE EQUATION 

To solve the Liouville equation to all orders in the 
probe field amplitudes E and B we first transform to 
a rotating frame and specialize to steady state solutions. 
This gives a set of 25 algebraic equations which we cast 
into a matrix form by arranging the 5 diagonal pn . . . P55 
and 20 off-diagonal density matrix elements p~2i ■ ■ ■ into 
a 25-dimcnsional vector p. We end up with an inhomo- 
geneous matrix equation 



fields E and B. At first glance the separation does not 
seem to be unique. To determine a EE (E, B), a EB (E, B), 
a BE (E, B), and a BB (E, B) numerically we formally ex- 
pand / and g in a power series in E and B 



f(E,B)=Y^f nm E n B r 



n.m 



9(E, B)=22 9nmE n B r 



(A.5) 
(A.6) 



To separate the electric and magnetic properties uniquely 
we use the fact that physically there must be an odd 
power of field amplitudes. In fact all but one (fastly rotat- 
ing) factors e~ luJpt must be compensated by factors e luJpt . 
Otherwise the (untransformed) polarizabilities would not 
oscillate with the probe field frequency u) p . Since an odd 
power of field amplitudes can only be realized by an odd 
power in E and an even power in B or vice versa an even 
power in E and an odd power in B wc formally split into 
even and odd subscries 



f(E, B)=J2 f* m \E\ 2n \B\ 2m E + £ f B JE\ 2n \Bf™B, 



Mp = a 



(A.l) 



1 {E,B) = Y,9nJn 2n \B\ 2m E + Y J 9n m \E\ 2n \B? m B. 



where the inhomogencity is given by the 25-dimcnsional 
vector a = (1,0,0,... ). The matrix M. contains all cou- 
plings, dctunings and decay rates for the system in ques- 
tion. To solve equation (|A.1[) for the sought density ma- 
trix elements ^34 and p~2\ we have to invert M. which can 
only be done numerically after specifying explicit num- 
bers for all parameters. 

In general we find that P34 and pi\ are functions of 
both, the electric and the magnetic field amplitude 



p 34 = f(E,B) 



p 21 =g(E,B). 



(A.2) 



We emphasize that the analytical form of the functions 
f{E, B) and g(E, B) is unknown. As we want to compare 
with the result of linear response theory we have to bring 
(|AT2|) in the form of Eq. OT 



(A.3) 



d 34 p 3i = a EE {E, B)E + a EB (E, B)B 



The polarizabilities are therefore given by the appropri- 
ate subseries, e.g. 

a EE (E,B) = d a4 Y,fnm\E\ 2n \Bf m , 



and similarly for a EB (E, B), a BE (E, B) and a BB (E, B). 
The utilization of symmetry properties then gives 



a EE (E, B) = ^[f(E, B) + f(E, -B)} (A.7) 



a EB (E,B) = ^±{f( E ,B) + f(-E,B)] (A. 



J (E,B) = Z±\g( E ,B)+g(E,-B)] (A.9) 



M21P21 = a BE (E, B)E + a BB (E, B)B 



BB 1 



(A.4) 



In contrast to the linear response theory here we deal 
with the exact solution of the Liouville equation and 
therefore the polarizabilities are still functions of the 



a BB (E, B) = ^[9(E, B) + g(-E, B)\ (A.10) 

which represents a unique numerical solution for the po- 
larizability coefficients. 
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